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Overview of Topics 

Pharmacometrics end Trtal Simulation 

Important Aspects 
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Population Based Modeling and 
Clinical Trial Simulation 

Chm.'ci Trial Simulation and Population Modeling are 

C'.ar.c iy arlaterf 

Trial Simulation Used to Qualify Population Based Models 

' vuvji fr. d .*t).e Check S.nvatc Data Using and 

mutated Data w-rh Observed Oat-j 

Components of Population Models and Simulation Models 
are Similar 
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■ - /'v.*c' Misti i.\y .fT-.u-c’ st t-t Can Be Uctd 
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Potential Uses of Both Must Be Defined end Specified 
Prior to Initiating Work 

•Vc:» * ■..cu.ifltii.iTs Moddc t jyw Used fer .Simulation 
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General Issues 

■: ;\:c £>;•.<;!ojimcrt Must Be Complete 

Covariate Model Must Be Developed 

• • ■ : •" •••’■''• '■ c •• r jP cr. rtdos-sunpttcr c 

PK/Pb Relationship Must Be Established and Confirmed 

v t'r r»:H *Y Pk'.-Pf* 

Outcome Model Must Be Built 

Dose Regimen Options Must Be Reasonable and Consistent 
with Clinical Practice 


Analysis Plan Must Be Developed Prior to Initiation 


Sufficient Resources Must be Given to Modeling and 
Simulation 


General Issues 

Hv/iv Should Drug Development Programs be 
fc-rs-orii'd to generate Sufficient Information About 
Advr- M sc Events and Outcome? 

Dose and Schedule Dependence 

(jp Vit or f, - v,nt 5nr.il ir atjrfitt 

for Surrogate Variables / Informative 

Co variate;, 

Genetics or Genomics Data 
Learning and Confirming 

«*■':** ^5 th<* Forty Pho'.n ror be Used or.d Tested 

Noe{l tc Collect Mere ond Better Information on 
/>dv.'-"“e G'.ort: and Response 

Atificial Inteilegence and Fuzzy L091C Methods 

vr-rcuj/ .■jh-'.iv' rc- H-ive Util ty to ftecc§r.,p.:- 3rd ^dvcric 

Collection of Patient Diorv Information or Other Means 


Selection of Desired Product 
Using Cost Functions 

Surface Methods or Combinatorial 

>p*-?lVZ AliW 

Describe the Relationship Between Desired Product I 

Performance and Costs 

’ yr■ - '• :>• o ■.">> j ;P r- ft *i orbing Benefits ogoinst Cosis. eg. Less 

L •v:}. , nt 3 Ajj T..t Increased Crug c’o-t 


Uses Objective Function To Select Optimal Solution 
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Basic Questions': f 

Basic Questions to /Ask Prior to Beginning a 
Modeling or Simulation Analysis 
■ What bo You Need to Know? 

The design of a simulation study is as critical as tKe 
design of a clinical protocol. \ 

How Well Do You Need to Know it? J 

Many models must be defined for simulation 
Need to determine which ones must be "correct" arjd 
which ones can be “OK" but not perfect 
What Are the Consequences if You Are Wrong? 
Hypothesis testing or sensitivity analysis must be ! 
conducted. 

Models must be qualified 
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Population Based Modeling and Sirnutotipn 
Terminology 

Need to Define Several Models 

Covariate Distribution Model 

Defines the characteristics of the patients frr the simulation experiment I 
« Demographics 

- Concomitant disease or disease stage 

Input / Output Model 

Defines the clinical indication and the drucf that is b</nff tested 

■ Structural and Stochastic Models 

* PK, PD, DJkom Prapfesslen, Outcomes 

■ ftsvonate Models 

» Defines effect of ccyorlctes tn PK, PD g* Outcome Structure Models 
Execution Model 

Everything Else! ' 

■ Dosing Practice awt Comjtfanc* 

» Drop out 

• Study Evaluation 


Mode! Components " 



Age 1 

A C c2 

^WeiRfct 2 
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Covariate Distribution Models 
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Covaria+e 

Joint Function 

Joint Function Describes | 
the Relationship of 
Patient CflvariateS (e.g. 

Age and CrCL) 

Establishes a 
Multivariate Density 
Function of These 
Covarrates 

"Virtual Patients" Will 
have Reasonable 
Covariate Combinations 


Other Methods for Dealing with 
Missing Data 

Case Deletion ” ■ 

‘ Delete Patients With Missing Data 

If demographic data, remove entire patient \ 
If PK observations, delete missing data items 
Problem! If Many Patients Have Missing Data,' 
Data Interpretation Can be Incorrect 

Covariate or demographic model information: may be; 
lost 

If many observations are below the limit of assay . 
quantification, then modeling may be difficult 
If patients dropped out, assessment of outcome is 1 
difficult 


Other Methods for Dealing wifR'TT ? I 
Missing Data .; _ ..r'Tvtl 

Imputation 

• Replace Missing. Demographic Values With Median Valuq. 

Determined From the Rest of the Data ■ 

Used for missing covcriatcj - easy to do if 

Problem! although thte Imputation will maintain the median vatue 
of the distribution, it niters the secondary characteristics of the 
distribution {*. 3 - variance^ 

Replace Missing Observations Due to Assdy Limits with ^ 

Lower Limit of Assay Quantitation 

This is a fairly robust technique if data are relatively rich ; 

Fit the Data » 

Only useful if the number of missing observations is low I 
(approximately Jess than 5% of the total' number of observations) 1 
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Case Deletion 

■- Easy to do but can alter data interpretation '' 

Imputation 

• Easy to do but limited utility and can after data 

interpretation ” ;f" 

Multiple Imputation 

• Difficult to do but can be useful, particularly 
useful for non-measurable PK or PD values 

Joint Function j 

• Somewhat difficult to do but most robust methjad 

for replacing, missing data especially for ; 

covariates , 
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Pharmacokinetic (dose, concentration, time) ' 

* drug disposition in individuals A populations ! 

■ disease state effects (renal A hepatic dysfunction) 
intervention effects (hemodialysis) 
concurrent medication effects 
pharmacogenetic influences 

Pharmacodynamic {dose or concentration, effect, 
time) 

physiologic A bicmerkers 

• surrogate endpoints 

■ clinical effects and endpoints 
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Empirical and Mechanistic PK/PD 
Models 

(Sole of Empirical PK/PD Models 

■ Sometimes Empirical Approach is the Only Option * 

• Specialized Applications 

Results May Guide Selection of Mechanistic Model 

■ Computationally Easter to Develop 
Role of Mechanistic PK/Pb Models 

Prior System Information Useful ‘ 

Knowledge Gained on Drug Effectcs has General Utility 
Extrapolation and Prediction of New Regimens or ? 
Formulations More Reliable 
■ Numerically Difficult To Develop and Parameter 
Identif iability Can be Problematic 
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Disease Progress Model Building Process 

Talk to a Disease Specialist } 

Draw pictures of time course of disease.! 

Translate into disease progress model ! 

Explain the models/parameters I 

Ask Disease Specialist for advice on 
factors influencing parameters I 

Translate into models of parameters ) 


Example Construction of a Disease 1 
Model 


Cadaveric. Dow; 

Mulched IK* Uo mulched? 
ttrsi Transplant? 



AtlrninjiUT 
Druffor Placebo 



Measure CD25+T Cells Measure 11.6, TNFalpba 


Data for Disease Progress Models” 

Use Measured Concentrations and Observations r 
• This Is easy To do ' . "C 

Use a “Link" model or Effect compartment) to 
Create Lag Between Observed Concentrations and 
Effects r 4 : 

This can be helpful os the time course for change in ! 
disease status is usually not the same as the time course 
of the drug 

Use Probability Models if no Continuous 
Measurement Available 
Literature 
■ Expert Input 


PM3001070320 


Source: https://www.industrydocuments.ucsf.edu/docs/rrjk0001 



















Issues Involving Building Disease 
Progress Models 

Lack of Continuous Measurements for 
Building Disease Progress Model 
» Literature Data 

Expert Input ; 

Resources 

Lack of Available Data for Untreated Patients 
■ Time Required to Collect Data and Develop Model 
• Variability Inherent in Data May Require Large 
Numbers of Subjects to Determine Parameters 
Accurately 


Role of Modeling Clinical ! 
Outcomes 

Improve Mechanistic Understanding . i 
Detect Factors Predictive of Inter- and 
Intra- Patient Variability for Both PK 
and Adverse Events 

Predict Consequence of Changing PK, PD 

and Dose Regimen 

Dose Schedule Optimization 


Clinical Outcome Models 


i .«I 


Tmituwiit OuKopw 
P ativiil t 
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Adverse Events 
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Adverse Event Classification 

Adverse Event Relationship to Test Drug 

Classifications; : ' 

Definitely. Probably, Possibly, Unlikely 
trades . ... 

None, Mild, Minor. Moderate, Severe 
Frequency and Immediacy of Adverse Event 
"A" Classification (Accountable, Predictable) 

• High Frequency, Immediate 
High Frequency. Delayed 
Low Frequency, Immediate 
Low Frequency. Delayed 

• "B ,r Classification (Bizarre, Unpredictable) 
Idiosyncratic 
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PK/PD Relationships for iHinical 
Outcomes 


.s«iniwr> PK vuiijblc (<uil L'Krlinc-coucw 


luU On i come font* 
IllU&Crtlifte: 


Highest Grade 

Highuit Crnuic Versos Cp-t 

Versus AUC 


Outcome 

Outcome 

Time Course 

Time Course 

Versus AllC 

Versus Cp-t 
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Why Use the Full Time Course' of 

The Clinical Outcome Information? 

Maintain All Information. 

More Than a Single Aspect of the Outcome 
May be Important, ; 

Time Course of Outcome May be Correlated with 
Concentration : 

Individuals with bifferent Quantity/Quality 
of Data Can be Pooled. 

Population Based Approach 
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Timing of PK/Outcome Mocfel j" 

Development ) 

Continuous PK and Continuous Outcome Data 
■ PK/AE Models Can be Started in Early Phase of 
Development I 

Continuous PK and Summary Outcome Data 

Usually Require Large Numbers of Subjects to Avoid Bias 
Often Need to Wait Until Phase XI j 

Pooling Data Across Indications May be an. Option ' 

Summary PK and Summary Event Data j 

Usually Require Large Numbers of Subjects to Avoid Bias 
Often Need to Wait Until Phase n 
Pooling Data Across Indications May be an Option 


Integrating PK/Outcome Models , " 
into Drug Development ' 

Development of PK/Outcome Model is Frequently 
Conducted Retrospectively " 

■ Integrated PK/AE and PK/Efficacy Analysis'I’rovIdes 
Opportunity to Test Model and Improve Drug Dose 

Regimen . ' 

Approximately 20% of Marketed Drugs - Undergo Post- 
Marketing Dose Adjustment! 

Integrated approach Lends Statistical and Scientific 
Validity to PK/Outcome Model Results 
Use of PK/Outcome Models Can Provide Better ‘ 
Assessment of Treatment Effects Than Traditional 
Statistical Assessments 

Permutation or Randomization Test 
Utility for Clinical Trial Simulation 


Approaches for Analyzing Summary*^ 

PK/Clinical Outcome Data " ; T" 

Approach Dependent on Data Available and 
Project Needs ' . " , 'Sif t; 

■■■ Continuous PK/Outcome Model 1 

! Categorical or Ordered Categorical Logistic '' 

Regression 

Adverse Event or Efficacy Data as Available only as 1 
Srades i ■ 

- Hazard Models 1 

Assessment of Risk of event or survival. 

• Time to Event Models 
Count Models 

Assessment of frequency 
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Development of Continuous P£7butco”me 
Models ! 

Determine the Type of PD Model That Best Describes 
Response 

Direct Effect j". “ 

No Lag Time to 8ff*et 

Stimulatory or Inhibitory Effect? ' 

p«sejir>e Response? \ 

Effect Compartment 

Short Leg Time due to Diffusion? 

Stimulatory or Inhibitory Effect? 

Baseline Response? 

• Indirect Effect? 

Longer Lag Time Dae to Drug Effect Occurring Many Steps Prior to 
Event 

Stimulatory or Inhibitory? 


Diagnostic Plots for Continuous Evelrff 
Data 


Plot of Concentration vs Effect No hysteresis Is evident 
so a dlrecteffecjt model tan be used. This drug appears 
to be stimulatory !..' . 


Concentration 


Plot of Concentration vs Effect. Hysteresis isevicfcnt 
so an effect compartment model can be used. Thfe 
drug appears to be stimulatory . i 


Concentration 


$• •- 




Time 


Plot of Concentration a nd Effect vers us time, There 
Is no effect until drug. Is goo®from the patient. A. 
hysteresis plot would not work here because the lag 
time is so long. Anindlrecteffect model should be used; 
This drug appears to be stimulatory 


Paclitaxel Pharmacokinetics 


Non-tinear elimination A non-linear 
distribution (Sonnichse n et al 1994) 

i Tj(T.^n 

ED 


Periplwui 



ED 
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% cj incline 


Paclitaxe! Neutropenia 





Cell Model for Hematologic Toxicity 


FectUwri; 

■ "' TOte 

CftU '• 


: Kilcji 


l-rflia: Fi'IbCt's, )krULlttk*<on,»Mi, SwuWrtint.Spyrrrtwwm, Knrtewiu l« prcpiuatka. - 


PM3001070327 


Source: https://www.industrydocuments.ucsf.edu/docs/rrjk0001 













Predicted Schedule Dependence ' " 



Threshold Models of /Adverse Event? 


jn J» ii> 
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Comparison of Threshold Model with full 
Model f 



Categorical Data 




Measurement Scale is a Set of Categories 1 ’’ 

■ Nominal (c.g, Success / Failure) ’ ^ 

■ Cardinal (e.g. Caucasian Asian Hispanic) j . 

Non Ordered ond Unrelated !. 

• Ordinal (e.g. Mild Moderate Severe) 

More Than 2 Response Levels < 

Response levels are Related and Ordered r \' 

Mare Than One Observation Per Patent 
Usually Considered Qualitative But May be Treated 
as Quantitative 

m Levels and (m- 1 ) Indicator Variables 
Levels Used as Actual Values 


Logit Models for Categorical Data J ^ Jj 
Analysis of Ordered Categorical Response fs 
Usually Referred to as "Logistic Regression" 

For a Probability (re) Between 0 and I 
> Logit Model : g(jt) = log [it / (1-it) ] [ 

= a +■ [t 

Maps the Probability onto a + Infinity to - 1 
Infinity scale t 

Solving for the Probability : 
ti(X) = exp (a ♦ j3X) / (1 + exp (a ■* px ) 
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Building Logit Models 


i ■ 


Frequency Trends Must be Established 
' Plot of Adverse Event Frequency by Grade and. Treatment 
Models Sensitive to Initial Estimates 

■ Transform Raw Data to Logit Plots 

■ Use Estimates: from Graphic Evaluation , 

Generate Data Base of Summary PtC Data and 

Adverse Events 

Can Do Simultaneous Continuous / Odd Type Fit [f Desired 

Fit Data Using Laplacian with Likelihood Option | 

Testing: 

1 Plot of Raw Data versus Predicted Probability of Response 
Bootstrap to qualify model 


Estimated Exposures 
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Percentage of Patients with Incidence of Bfeeding vs. 

. Groups of Ascending AUC Values 
(n=94 for Group 1 (Placebo), n=45 for Groups 2-B and 
n s 35 for Group 9) 
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Continuous Pharmacokinetics and 
Categorical Clinical Outcome Data 

Options To Handle Data -P 

Simultaneous Fit of Continuous and Categorical 
Summary Data 

■ Serial Fit of Continuous and Categorical Summary 
Data • ' "• ' 

Investigate Various Summary Parameters for PIC for 
Predictive Ability 

Css, AUC, Cmax ! 

Selection of Analytical Approach Depends oni 

Needs of Development Project and on Available 

Data \ 


Hazard Functions 

Hazard Functions Define the Rate of 

Occurrence of An Event 

Instantaneous Progression 
> Hazard Can Be Defined Using Several Functions 
Constant 
Linear 
Exponential 
Other (/,*., Weibull) 


tf.r-tf I 


: THP 
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Hazard Functions 

Cumulative Hazard is the Integral of the 
Hazard Over a Pre-Defined Period of Time 
Describes the Risk 

■ Translates or Correlates Pharmacodynamic 
Response into a Useful Measure of Outcome 
Assessment of Likely Benefit or Adverse Event 
Comparison With Existing Therapy 
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Survival Function 


'"T 


S(t) = P(T > t) 

Monotone, Decreasing Function 
Survival is 1 at Time=0 and O as Time 
Approaches Infinity. 

The Rate of Decline Varies According to Risk 
of Experiencing an Event | 

Survival is Defined as i 

S(t) ■- exp(-/y(/» 


Using Hazard FunctiorisTih 

PK/PD Models 1! ^ 

If Hazard Function is Defined as a Constant " 
Rate "K" Such that r T ‘ 

h(t) =K ' \ - 

Then the Cumulative Hazard is . '"V 

H{t\ - \Kdt '= Kt I 


Kt =~ in [££;}]: 


Survival is 


S(t) = exp [- Kt] 


Hazard, Cumulative Hazard, 
and Survival 


In This 

Example 

Hazard 

Remains 

Constant 

Cumulative | 

Hazard 

Increases 

With Time} 

Survival 

Drops 


ml 
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Example: Comparing Neutropenia With 
Different Hematopoietic Factors 



Comparing Neutropenia with Differghtl 
Hematopoietic Factors 
New Product Reduces Nadir and Recovery " 
Time ' 

1 Is the Improvement Enough to See Clinical 
Benefit? 

Evaluation of New Product More Relevant if 
PK and Pb Changes Related to Clinical 
Outcome 
■ Hazard Functions 


Comparing Hematopoietic Factors Us ing 

Hr ■“ - - 















